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Abstract 

We present two new analytic formulations of the Density Matrix Renormalization Group Method. 
In these formulations we combine the block renormalization group (BRG) procedure with Variational 
and Fokker-Planck methods. The BRG method is used to reduce the lattice size while the latter 
are used to construct approximate target states to compute the block density matrix. We apply our 
DMRG methods to the Ising Model in a transverse field (ITF model) and compute several of its 
critical properties which are then compared with the old BRG results. 
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1 Introduction 



The Density Matrix Renormalization Group Method (DMRG) introduced by White fl]] has arguably 
become the most powerful numerical tool to retrieve the essential features of interacting quantum 
Hamiltonians in ID such as spin systems [Q], Q and fermion systems Q. This method is a real-space 
renormalization group (RG) method which is specially well suited when dealing with zero temperature 
properties of many-body systems, a situation where the Quantum MonteCarlo methods happen to be 
particularly badly behaved as far as fermionic systems are concerned Q. 

The origin of the density matrix RG method relies on the special treatment carried out by White 
and Noack || on the ID tight-binding model, the lattice version of a single particle in a box. It was 
Wilson the first to point out the relevance of this simple model in understanding the sometimes bad 
numerical performance of the standard Block Renormalization Group (BRG) method. In reference || 
the authors proposed a method called Combination of Boundary Conditions (CBC) which performs 
extremely well as compared to the exact known solution of the model. Recently, we have clarified 
the role played by the boundary conditions in the real-space renormalization group method || by 
constructing a new analytical BRG-method which is able to give the exact ground state of the model 
and the correct l/iV 2 -law for the energy of the first excited state in the large A^(size)-limit. 
The problem with the CBC method is the difficulty when trying to generalize it to interacting models. 
In H yet another method was presented called the Superblock Method which is the precursor of the 
density matrix RG method. The DMRG method is able to treat many-body problems as for example 
the ID Heisenberg model of spin S = 1 |ffl|, ^ which happens to be a non-integrable model. 
In this paper we present two novel versions of the DMRG method based upon the Perturbative- 
Variational and Fokker-Planck approaches to quantum lattice Hamiltonians recently introduced in 
references |7| and [Q respectively. We arrived at these new methods by searching for an analytical 
formulation of the density matrix RG method. As it happens, the usual implementations of the DMRG 
method are intrinsically numerical for they rely on Wilson's procedure of enlarging the system size in 
his RG-treatment of the Kondo problem Q . We find interesting to address the problem of constructing 
analytical extensions of the DMRG method for several reasons. Firstly, it is known that standard block 
RG methods proved to be useful when dealing with qualitative features of some important models such 
as the ITF model [{Tcfl , Lattice Gauge Model [|ll]], Heisenberg model [12|, Hubbard model jj~3[l , etc. 
and we want to see how the DMRG method performs when compared to those analytical BRG treat- 
ments. To this purpose, we have applied our Variational and Fokker-Planck DMRG methods to the 
Ising model in a tranverse field (ITF model). There are more interesting models, but the ITF model is 
simple enough for a first application of these methods. As a matter of fact, it was Drell and the SLAC 
group [[H]] who started to apply the standard Block Renormalization Group menthod to study QCD 
and they used the ITF model as a test model. Later on in the 80's this BRG method has been also 
applied to the study of strongly correlated systems and its implications in High-T c Superconductivity. 
Secondly, the DMRG method as it stands is also intrinsically one-dimensional and it is an open prob- 
lem to find feasible numerical schemes to work with higher dimensional systems. As it happens, our 
Variational and Fokker-Planck DMRG methods are generalizable to dimension higher than one and we 
might consider them as a first attempt to solving this important extension of the DMRG method. 
This paper is organized as follows. In Sect. 2 we present a brief introduction to the Block Renormal- 
ization Group methods based upon the concept of the intertwiner operator T. This allow us to make 
a unified formulation of both the standard BRG methods and the new Density Matrix RG method 
according to our analytic formulation. In Sect. 3 the Block RG-method is applied to the ITF model 
following the formulation of Sect. 2 and several critical exponents are computed. In Sect. 4 we present 
our new DMRG methods and apply them to compute the intertwiner operator and critical exponents 
for the ITF model. The results are compared with the old BRG results and we find that the density 
matrix methods perform better. Sect. 5 is devoted to conclusions and prospectives. 



1 



2 Block Renormalization Group Methods (BRG): Brief Review 



In this section the block renormalization group method is revisited and we present a new and unified 
reformulation of it based on the idea of the intertwiner operator T to be discussed below. This for- 
mulation will allow us to introduce the new Variational and Fokker-Planck DMRG methods on equal 



footing as the standard BRG method. For a more extensive account on this method we refer to [14] 
and chapter 11 of reference and references therein. 

The block RG-method is a real-space RG-method introduced and developed by the SLAC group 
[10]. Let us recall that Wilson developed his numerical real-space renormalization group procedure 



to solve the Kondo problem |9J. It was clear from the beginning that one could not hope to achieve 
the accuracy Wilson obtained for the Kondo problem when dealing with more complicated many-body 
quantum Hamiltonians as the ones mentioned in the introduction. The key difference is that in the 
Kondo model there exists a recursion relation for Hamiltonians at each step of the RG-elimination of 
degrees of freedom. Squematically, 

Hn+i = Hn + hopping boundary term (2-1) 



H N+1 = R(H N ) (2.2) 

The existence of such recursion relation facilitates enormously the work, but as it happens it is specific 
of impurity problems. 

From the numerical point of view, the Block Renormalization Group procedure proved to be not fully 
reliable in the past particularly in comparison with other numerical approaches, such as the Quantum 
MonteCarlo method which were being developed at the same time. This was one of the reasons why 
the BRG methods remained undeveloped during the '80's until the begining of the '90's when they are 
making a comeback as one of the most powerful numerical tools when dealing with zero temperature 
properties of many-body systems. 

Let us first summarize the main features of the real-space RG. The problem that one faces generically 
is that of diagonalizing a quantum lattice Hamiltonian H, i.e., 



H\4>>=E\iP> (2.3) 

where \ip > is a state in the Hilbert space TC. If the lattice has N sites and there are k possible states 
per site then the dimension of 7i is simply 

dimH = k N (2.4) 

As a matter of illustration we cite the following examples: k = 4 (Hubbard model), k = 3 (t-J model), 
k = 2 (Heisenberg model) etc. 

When N is large enough the eigenvalue problem (2.3) is out of the capability of any human or computer 



means unless the model turns out to be integrable which only happens in some instances in d = 1. 
These facts open the door to a variety of approximate methods among which the RG-approach, specially 
when combined with other techniques (e.g. numerical, variational etc.), is one of the most relevant. 
The main idea of the RG-method is the mode elimination or thinning of the degrees of freedom followed 
by an iteration which reduces the number of variables step by step until a more manageable situation is 
reached. These intuitive ideas give rise to a well defined mathematical description of the RG-approach 
to the low lying spectrum of quantum lattice hamiltonians. 

To carry out the RG-program it will be useful to introduce the following objects: 



• H : Hilbert space of the original problem. 
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• Ti': Hilbert space of the effective degrees of freedom. 



• H: Hamiltonian acting in TC. 

• H': Hamiltonian acting in TC' (effective Hamiltonian). 

• T : embedding operator : TC' — ► TC 

• :truncation operator : TC — > TC 

The problem now is to relate H, H' and T. The criterium to accomplish this task is that H and H' 
have in common their low lying spectrum. An exact implementation of this is given by the following 
equation: 

HT = TH' (2.5) 

which imply that if *$>' E , is an eigenstate of H' then T*$>' E , is an eigenstate of H with the same eigenvalue 
(unless it belongs to the kernel of T: T^>' E , = 0), indeed, 

HT^' E , = TH'ty E , = E'T^' E , (2.6) 
To avoid the possibility that T^' = with \P' ^ 0, we shall impose on T the condition, 

rtr = i w (2.7) 

such that 

* = IV ^ itf = rt^r (2.8) 



Condition (2^7) thus stablishes a one to one relation between TC' and Im(T) in TC. 



Observe that Eq. ( |2.5D is nothing but the commutativity of the following diagram: 

TC' ^ TC 
H' I I H 

TC' TC 

Eqs. (|2.5| ) and (|2.7| ) characterize what may be called exact renormalization group method (ERG) 
in the sense that the whole spectrum of H' is mapped onto a part (usually the bottom part) of the 
spectrum of H. In practical cases though the exact solution of Eqs. fl2.5| ) and (|2.7| ) is not possible so 
that one has to resort to approximations (see later on). Considering Eqs. Q2.5[ ) and ( p. 7] ) we can set up 
the effective Hamiltonian H 1 as: 

H' = T^HT (2.9) 

This equation does not imply that the eigenvectors of H' are mapped onto eigenvectors of H. Notice 
that Eq.( |2.9| ) together with ( |2.7| ) does not imply Eq. (|2.5|). This happens because the converse of 
Eq.(|2jD, namely TT* + l H 

is not true, since otherwise this equation together with (|2.7| ) would imply 



that the Hilbert spaces H and H' are isomorphic while on the other hand the truncation inherent to 
the RG method assumes that dim H' j dim H. 

What Eq.( |2.9[ ) really implies is that the mean energy of H' for the states of TC' coincides with 
the mean energy of H for those states of TC obtained through the embedding T, namely, 

< ^'\H'\^' >=< TV'\H\TV' > (2.10) 
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In other words T^' is used as a variational state for the eigenstates of the Hamiltonian H. In 
particular T should be chosen in such a way that the states truncated in TC , which go down to H! , are 
the ones expected to contribute the most to the ground state of H. Thus Eq. ( p.9Q is the basis of the 
so called variational renormalization group method (VRG) []. As a matter of fact, the VRG method 
was the first one to be proposed. The ERG came afterwards as a perturbative extension of the former 
(see later on). 

More generally, any operator O acting in 7i can be "pushed down" or renormalized to a new operator 
O' which acts in TL' defined by the formula, 

O' = T ] OT (2.11) 

Notice that Eq.( |2.9| ) is a particular case of this equation if choose O to be the Hamiltonian H. 

In so far we have not made use of the all important concept of the block, but a practical implemen- 
tation of the VRG or ERG methods does require it. The central role played by this concept makes all 
the real-space RG-methods to be block methods. 

Once we have established the main features of the RG-program, there is quite freedom to implement 
specifically these fundamentals. We may classify this freedom in two aspects: 

• The choice of how to reduce the size of the lattice. 

• The choice of how many states to be retained in the truncation procedure. 

We shall address the first aspect now. There are mainly two procedures to reduce the size of the lattice: 

• by dividing the lattice into blocks with n s sites each. This is the blocking method introduced by 
Kadanoff to treat spin lattice systems. 

• by retrieving site by site of the lattice at each step of the RG-program. This is the procedure 
used by Wilson in his RG-treatment of the Kondo problem. This method is clearly more suitable 
when the lattice is one-dimensional. 



We shall be dealing with the Kadanoff block methods mainly because they are well suited to 
perform analytical computations and because they are conceptually easy to be extended to higher 
dimensions. On the contrary, the DMRG method introduced by White []l|] works with the Wilsonian 
numerical RG-procedure what makes it intrinsically one-dimensional and difficult to be generalized to 
more dimensions. Thus why we shall formulate our Variational and Fokker-Planck DMRG procedures 
as block renormalization methods in section 4. 

Block RG-methods have recently received also renewed attention in one-dimensional problems in connec- 
tion to what is called a quantum group symmetry [16]. Based upon this symmetry we have constructed 



a new BRG-method that we call g-RG which among other features it is able to predict the exact line 
of critical XXZ models in the Anisotropic Heisenberg model, unlike the standard BRG-method. 

To exemplify the standard BRG-method we shall study a ld-lattice Hamiltonian, the Ising model in a 
transverse field (ITF model) . The main ideas are also valid in higher dimensions although computations 
are more involved. Hence we shall be dealing with a one-dimensional lattice, usually a periodic chain. 
In every site of the chain there are k degrees of freedom, hence: 

H = C k ® X. ®C k := ® N C k (2.12) 

We shall consider Hamiltonians H containing operators which involve only a single-site part H$ or 
two- nearest-neighbour-site part H ss and will be simbolically depicted as in Fig.l, in such a way that, 

1 The word variational here is used with a different meaning as in the introduction of this paper where it refers to the 
variational choice of the target state in the DMRG method to be discussed in section 4. 
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H = H S + H ss (2.13) 

As a matter of illustration, let us give one example of this decomposition in the ITF model which will 
turn out to be very useful in putting many key ideas to the test. 

3 Block RG- Approach to the Ising Model in a Transverse Field (ITF) 

The Ising Model in a Transverse Field is originally a one-dimensional quantum lattice system with 
quantum critical properties equal to the well-known thermal critical properties of the classical 2D-Ising 
Model. The lattice Hamiltonian of the ITF model is: 

N N 

H N (T, J) = -TJ2 of - J £ o)o)+ x (3.14) 

3=1 3=1 

where a? and <rj are the standard Pauli matrices acting at the j-th site of the chain. 

The Hilbert space of states and the intrablock and interblock Hamiltonians for this model are, 
respectively: 

H = <C 2 (3.15) 

N 

H s = -Tj2°j (3-16) 

3=1 
N 

H ss = -jJ2° Z j Vj+i (3-17) 

3=1 

The first step of the BRG method consists in asembling the set of lattice points into diconnected 
blocks of ub sites each, as in Fig. 2. 

In this fashion there are a total of N' = N/ub blocks in the whole chain. This partition of the 



lattice into blocks induces a decomposition of the Hamiltonian (2.13) into an intrablock Hamiltonian 
Hb and a interblock Hamiltonian Hbb as illustrated in Fig. 3. 

Observe that the block Hamiltonian Hb is a sum of commuting Hamiltonians each acting on every 
block. The diagonalization of Hb can thus be achieved for small ub either analytically or numerically. 
The content of Fig. 3 can be written as 

H = Hb + XH BB (3.18) 

where A is a coupling constant which is already present in H or else it can be introduced as a parameter 
characterizing the interblock coupling and in this latter case one can set it to one at the end of the 
discussion. 



Eq. ( 3,1§| ) suggests that we should search for solutions of the intertwiner equation (^j) in the form of 



a perturbative expansion in the interblock coupling constant parameter A, namely, 

T = T + ATi + A 2 T 2 + . . . (3.19) 

H' = H' + \H[ + X 2 H' 2 + ... (3.20) 
To zeroth order in A Eq. ( |2.5D becomes 

H B T = T H' (3.21) 
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Since H B is a sum of disconnected block Hamiltonians h -, , f = 1, . . . , N implicitly denned through 
the relation 

N' 

H B=Y, h ? ( 3 - 22 ) 
i'=i 

one can search for a solution of To in a factorized form 

N' 

t o = n r o,/ ( 3 - 23 ) 

and an effective Hamiltonian H' which acts only at the site f of the new chain, 



'o 

N' 



H' =J2 h f ] = H 's' (3-24) 
j'=i 

Observe that H' s , is nothing but a site-Hamiltonian for the new chain. Eq. ( |3.21| ) becomes for each 
block: 



hf>T 0if = T 0if hf, > (3.25) 



k' h ns -k' 



The diagonalization of h-, for j' = 1,...,N' will allow us to write 



i=l a=l 

where |i >,■' for j = 1. . . . , k are the A; -lowest energy states of h), . Moreover, we suppose that h), is 
the same Hamiltonian for each block so that does not depend on the block. 

(s) 

The truncated Hamiltonian h -, and the intertwiner operator Tq ~< are then given by: 

j 

W 

1=1 
fe' 

r o,i' = e j'(*r ( 3 - 2§ ) 

i=l 

Later on we shall show examples of these relations. 

To obtain the first order correction to the Hamiltonian H[ we must consider Eq. ( [2.5D to first order 
in A: 

H BB T + H B T X = T H[ + TW (3.29) 
Multiplying the left hand side by Tq and using TjTo = 1 along with H b Tq = TqH' we readly obtain: 

T^HbbTq + H^T! = H[+ TfcH'o (3.30) 

We would like to kill the term proportional to TqTl. For this purpose Eq. ( |2.7| ) which implies 
TqTi + TTq = 1 is not very useful. A resolution of this problem can be accomplished if instead of the 
operator T one uses another operator T satisfying the defining equations: 

HT = TH' (3.31) 
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-.t 



rjr = i w (3.32) 



Then Tq = Tq and TqT\ = in which case Eq. ( 3.30[ ) simply becomes: 



H[ = TqHbbTo = H' s , s , (3.33) 

We can summarize these results saying that up to first order in A, the effective Hamiltonian H' can 
be obtained using simply the zeroth order intertwiner operator To (see Fig. 4): 

H '(up to order A) = H 's' + = T ^ Hb + XH bb)Tq (3.34) 

This is precisely the prescription of Drell et al. [1C]. 

The second order correction to H' can be obtained again from Eqs. ( |3.31| )-( |3.32| ) and is given by 

H' 2 = tIHbbTi (3.35) 



There is a close parallelism between the perturbative solution of Eqs. ( ^5| ) or ( |3.3lQ - ( p. 32 ) and the 



pertubation theory of the Schrodinger equation for a Hamiltonian of the form Hq + \H\ . As a matter 



of fact, the normalization condition (3.32) for operators is equivalent to the standard normalization for 
wavefunctions < ^ol^A) >=< ^ol^o >= 1 that is adopted to avoid normalization complications. In 
what follows we shall mainly concentrate on the first order solution Eq. ( [3.34 ). 



The final outcome of this analysis is that the effective Hamiltonian H' has a similar structure to the 
one we started with, H. The operators involved in H', and H' s , s , may by all means differ from those 
of H s and H ss , but in some cases the only difference shows up as a change in the coupling constants. 
This is known as the renormalization of the bare coupling constants. When this is the case, one may 
easily iterate the RG-transformation and study the RG-flows. 

Let us summarize the RG-prescription we have introduced so far in Table 1. 
We have denoted this prescription by BRGl(n s ,/e') where n s and k' have been defined earlier and 1 
denotes that we are working to first order in perturbation theory. 

For the ITF model we shall consider blocks of two sites (n s = 2) and truncation to two states 
(k' = 2). The block Hamiltonian has two sites and has the form, 

= -T(af + al) - Ja{a z 2 (T, J > 0) (3.36) 
The eigenstates of this block Hamiltonian ( 3.36| ) are given in increasing order of energies by, 



1 



\G >= -==(|00 > +a|ll >) E = -v / J 2 + 4r 2 (3.37) 
VI + a 2 

\E >= 4=(|01 > +|10 >) E = -J (3.38) 



V2 V 
1 

7i v 



\E' >= -7=(|00 > -|11 >) E = J (3.39) 
1 



\E" >= -__(_o|01 > +|10 >) E = V J 2 + 4r 2 (3.40) 
vl + a 2 

|0 > and |1 > are the eigenstates of a x , 

a x \0>=\0> 1 a x \l>=-\l> (3.41) 
and a = a(g) is the following function of the ratio J/2r := g, 
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which in turn satisfies 



(3.42) 



o(0) = 0, a(oo) = 1 (3.43) 
The intertwiner operator whithin each block has the form 

T (a) = \G)(0\' + \E) <1|' (3.44) 

where |0 >' and |1 >' form a basis of states at each point of the new chain. The effective Hamiltonian 
H' up to order J can be computed from Eq. ( 3.34| ). 

Thus to get TC we have to study the renormalization of the various operators entering in its definition. 
Using djU) and (^37|)-( glc|) one obtains after some elementary algebra: 



T t a x T 



1-q 2 
2(1 + a 2 ) ( 

1 + a 



1 + 



V2(l + a 2 



(1 + a) 



■1 + 



1 - a 



2(1 + a 2 )" ' 2(1 + a 

rn\z Z rp _ (1 + a ) 2 _z _2 

The range of the indexes run as follows: 

j = 2f-l+p 



2\ a f 



(3.45) 
(3.46) 
(3.47) 
(3.48) 

(3.49) 



j = l, 



,N 



(3.50) 



/ = l,...,iV/2 p = 0,l 
Applying Eqs. (gi4^) - (gi48|) to the ITF Hamiltonian, one gets 



rt^(r, j)t = ae + h n/2 (V, f) 



where 



AE 



N 
~2~ 



1-a 2 



+ 



J (1 + a) 



(1 + a 2 ) 2 (1 + a 2 



r' 



(1 + a 2 ) 2(1 + a 2 ) 



(3.51) 
(3.52) 

(3.53) 
(3.54) 
(3.55) 



2(1 + a 2 ) 

The derivation of Eqs. ([3.45 )- ( 3.48j ) and ( 3.52j )- p?55| ) does not make use of Eq. ( |3.42 ) and hence 
have a more general validity. In other words, we can use the function a(g) as a variational function in 
order to construct better ground states in the spirit of the VRG. 
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Eq. Q3.42 ) is one of the numerous choices we can make. We shall consider later on other examples. 
Schematically we can set up the following relationship, 



RG-Prescription with n s = 2, k' = 2 J ~ ® . , . (3.56) 

I a(0) = 0, a(oo) = 1 

The physical properties of Hn(T, J) depend only upon the ratio g = J/2T. If < g < 1/2 one is 
in a disordered region characterized by a unique ground state with unbroken symmetry (< a z >= 0). 
li g > 1/2 the Zi symmetry associated to the operator Q = FJj &p which commutes with Hn for N 
even, is broken. This is the ordered phase which has two degenerate ground states corresponding to 
< a z >= ±m / 0. 

At g = gc = 1/2 the system is critical and belongs to the same universality class as the 2D- 
classical Ising model. The critical exponents can be defined in terms of the behaviour of the "quantum 
observables" as functions of g c — g- 

Most of the critical exponents can be computed from the properties of the RG-transformation. In 
the case of the ITF model the RG-transformation can be obtained from Eq,( 3.53j )-( |3.55| ): 



9 ■= — ■= m = ^ +a{ f )? 1 ^^ 2 (3-57) 

y 2F yy> 2l-a{gY-g{l-a{g)Y V ; 



For any function a(g) satisfying Eq. ( 3.56j ), this transformation has 3 fixed points g* = 0,g c ,oo. The 



fixed points g* = and oo are attractive and correspond to the disordered and ordered phases respec- 
tively. The fixed point at gc is repulsive and correspond to the critical point of the ITF Hamiltonian. 
The value of the function a(g) at g = g c is a function only of gc and does not depend on the particular 
prescription chosen, that is, it is a universal function : 



u( \+ ( ^ 2^T^2g^+2g c -l 

9c = R{9c) => a(sc) = ^—^ (3-58) 

3 + 2g c 

This equation implies in particular that whatever prescription is chosen, assuming that a is real, the 
critical value obtained from Eqs. ( |3 . 53| )- ( |3 . 55| ) will always be less than the exact value 1/2. 

9c < 5c XaCt = 1/2 (3.59) 

To get the value of gc for a given prescription one has simply to find the intersection of the function 
a(g) and the function, 

C 2^1^+2,9-1 n < o < 1/2 
f(g) = 1 3+2 9 J - 9 - 1 ' / (3.60) 

y ' \ g > 1/2 V ' 

as is shown in Fig. 7. 

The analysis of the RG-equations usually has to be done numerically. However, there is a great 
deal of information that can be retrieved without completely solving the RG-equations if we know 
wheather the successive coupling constants g n — ► g n +i = R{9n) increase or decrease during the 
iteration procedure. To this end it is convenient to introduce the the familiar beta function /3(g) of 



quantum field theory which in this context is [10], 



(3{g) := R(g) - g (3.61) 

In Fig. 8 we have plotted the beta function for the ITF model we are analyzing. A fixed point of the 
transformation occurs at values of g which reproduce themselves under the RG-iteration, i.e., they are 
the zeroes of the beta function: 

f3{g c ) = (3.62) 
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There are 3 fixed points for besides the two zeroes at g = and 0.39, g = oo is also a fixed point for it 
cannot be reduced by further iterations. 

There is additional qualitative information which can be extracted from the shape of the beta 
function /3(g). In particular, the sign of (3(g) is responsible for the stability character of the fixed point. 
When (3(g) < (> 0) this means that g decreases (increases) after one iteration and the resulting g' 
lies to the left (right) of the g we started with. 

The outcome of this RG-analysis can be summarized by saying that the fixed points = 0, oo are 
stable fixed points while g c = 0.39 is an unstable fixed point. 



Given the RG-transformation Eq. ( 3.57 ) we can compute several critical exponents and compare 



them with the exact results in order to check the accuracy of the method. 

Correlation Length Exponent v. It gives the behaviour of the correlation length in the vicinity 
of g c 

£ ~ (9 - 9c)- u (3.63) 

Under the RG-transformation ( |3,57j ) £ —* £' = £/2 which leads inmediately to an expression for v in 
terms of the derivative of R(g) evaluated at the critical point, 



1 In R'(g c 



v In 2 

From Eq. ( gjg ) we can evaluate R'(g c ) as a function of a c = a(g c ) and a' c = -#| c , 



(3.64) 



A T := R'(g c ) = 1 + 2g c ( l —^f + /f ^^ (3.65) 

l + a c (l + a c ) A 

In Table 2 we show the value of v obtained for different choices of the function a(g). 

Dynamical Exponent z. At the critical point where g' = g holds, the Hamiltonian changes by an 
overall factor which in turn defines the dynamical exponent z 



In order to get z we notice that 



Hence 



H C ^H' C = ±H C (3.66) 



1 V V 

¥ = ( 7 )c = ( r )c (3 ' 67) 



Kj) _ 1 , ln[(l + a c 2 )/(l + ^) 2 ] 



(3.68) 

In 2 In 2 V ' 



It follows from Eq. ( |3.68| ) and the positivity property of a c ( ^.56 ) that z is always less than the 
exact value, 

z < z exact = 1 (3.69) 
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Magnetic Exponent [3. This critical exponent is denned through the spontaneous magnetization 
M z in the ordered phase, 



M =< a) > (3.70) 
Above a c but close to the critical point we will have 

M~(g-g c f (3.71) 

Gap Exponent s. Using scaling arguments satisfied by the BRG method, the gap G behaves like 

G~{g- g c ) s (3.72) 



with 



s = vz (3.73) 



Equation (3.46) relates the magnetization M and the one obtained after the RG-transformation 



M = l t + a M' (3.74) 
v/2(l + ^) 



Combining Eqs.( 3.74| ), ( |3.71| ) and g' = R(g) we arrive at, 



ML = [R'( 9c )f = 2^= / + QC (3.75) 
" • /2(1 + a%) 



Using Eq,( |3.68 ) we are able to relate the critical exponents (3, v and z through the following scaling 
relation, 

P = \™ (3-76) 

This relation is valid for any choice of the function a(g) and therefore it is characteristic of using a 
block containing two sites. Observe that the exact exponents of the ITF model never satisfy this scaling 
relation ( 3,76] ). This shows the limitation of the block method when using a two-site block. 



4 Density Matrix RG Method: Analytic Formulation 

The Density Matrix RG-method (DMRG) is an improved version of the real-space renormalization group 
methods introduced by White [RJ as a further elaboration of the ideas concerning the Combination of 
Boundary Conditions method ||. 

The fundamental difficulty of the BRG method lies in choosing the eigenstates of the block Hamilto- 
nian Hb to be the states kept. Since Hb contains no connections to the rest of the lattice its eigenstates 
have inappropiate features at the block ends. The CBC method of White and Noack is a first attempt 
to solving this intrinsic problem. The rationale of this method was that quantum fluctuations in the 
rest of the system effectively apply a variety of boundary conditions to the block. The CBC method 
proved to be very effective for the simple single-particle problem studied by White and Noak, but it 
happens to be ill-suited to interacting systems. The importance of the CBC method relies more on the 
lessons we can learn from it rather than the specific technicalities pertaining the simple case where it is 
applied successfully. There are two main ideas in order to proceed towards density matrix RG-method, 
namely: 
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• The block is not isolated. 



• The eigenstates retained are not eigenstates of a unique block Hamiltonian Hb- 

The DMRG method is in a sense an evolution of the CBC method in which we "let the system" to 
choose the best boundary conditions. White suggests that for a system which is strongly coupled to the 
outside "universe" (the rest of the lattice) , it is much more appropiate to use eigenstates of the block 
density matrix to describe the system (block), rather than the eigenstates of the system's Hamiltonian 
Hb- White's proposal can be stated by saying: 

• Choose to keep the ub most probable eigenstates of the block density matrix. 

It is possible to show that keeping the most probable eigenstates of the density matrix gives the most 
accurate representation of the state of the system as a whole . This is the basis of the Density Matrix 
Renormalization Group (DMRG) method. 

For the sake of completeness we present a brief introduction to the DMRG method. A superblock 
is called to a large block which contains several smaller blocks, one of which is the block to be used in 
the blocking procedure of the BRG method. Let us suppose that we have diagonalized a superblock 
and thereby obtained one particular state which is called the target state and probably will be 
the ground state. Let {\i), i = 1, . . . , 1} be a complete set of states of the the block B which we call 
"the system". Let also {\j), j = 1, . . . , J be a complete set of states of the superblock which we call 
"the universe" (see Fig. 5). Now we proceed to decompose the target state \tp) into its system- and 
universe-parts according to the following equation, 

IV>)=£^>b') ( 4 - 77 ) 

Next we want to devise a procedure to produce a set of states of the system denoted by 

\u a ), a = l,...,n B with \u a )=Y J uf\i) (4.78) 

% 

which are optimal for representing the target state in a sense to be specified below. The number of 
states kept is such that ub < I so that \ip) is represented approximately, that is, 

|V>}^}:=£ a "> a > \J) (4.79) 

where a a j are components to be determined by demanding that the following distance to be a minimun: 

V:=U)-\^)\ 2 (4.80) 
This minimization problem requires to vary over all a a j and u a subject to the condition, 

(u a \u a ') = <W (4.81) 

White has proved that the solution to this minimization problem is given by the optimal states u a being 
eigenvectors of the reduced density matrix of the system as part of the universe whose eigenvalues are 
the ub largest in mag nitude Q. 

The reduced density matrix for the system depends on the state of the universe which in this case is a 
pure state \tp). Therefore, the density matrix for the system is given by, 

PiV = ^i>iji>i'j (4-82) 

3 
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To summarize, when the entire lattice is assumed to be in a pure state, the optimal states to be kept are 
then ns most significant states of the reduced density matrix of the block, say B. 

Next step in putting density matrix RG ideas at work is to devise an efficient algorithm based upon 
these developments. A density matrix algorithm is defined mainly by two features, the same as the 
standard BRG algorithm is, namely, according to the form of the superblock and the manner in which 
the block are enlarged, e.g., doubling the block B' = B B (Kadanoff) or adding a single site B' = B + 
site (Wilson in Kondo). 

As far as computer power is concerned, generally is more efficient to enlarge the block by adding a single 
site rather than doubling the block. The reason for this is that the diagonalization of a superblock 
composed of say p identical blocks is difficult for a many-particle interacting system for the dimension 
of the Hilbert space of states goes like n p B , assuming that ns states are kept per block. 
This has led White to propose a variety of density matrix algorithms for both finite and infinite size 
systems which rely on Wilson's method for enlarging the system. Thus, they are intrinsically one- 
dimensional methods. 

In this paper we are presenting two new DMRG methods (Variational and Fokker-Planck) which are 
based on the blocking procedure of reducing the system size as the standard Block Renormalization 
Group method ot the previous section. This makes them potentially well-suited to address higher 
dimensional systems. 

Moreover, the DMRG algorithms introduced by White are intrinsically numerical for they are based 
upon Wilson's procedure of enlarging the system size. On the contrary, the BRG study of the ITF 
model carried out in the previous section was done in an analytical fashion mainly because it is a 
blocking procedure. We shall present hereby a new treatment of the ITF model along the lines of the 
density matrix RG method which is analytical. 

To this purpose, we have incorporated the block method in the DMRG algorithm. In addition, we shall 
incorporate another ingredient for the algorithm to become analytical: we shall choose the target state 
(ground state) using a variational method or a Fokker-Planck method, and then compare the results 
obtained in these two fashions with the standard BRG results of the previous section. 

The way to combine the above tools is as follows. We first set up a variational ground state for the 
whole chain whose energy is determined by solving the corresponding minimization equations. Next 
we use this state to construct a block density matrix p for a block having two sites in the philosophy 
of the BRG method. This p turns out to be (see below) a 4 x 4 matrix whose two largests eigenvalues 
denoted by \A) and \B) are kept to construct the intertwiner operator To as in sections 2 and 3. Once 
we make contact with the blocking method, the iteration procedure goes over and over. The important 
point now is that we can keep the DMRG study at an analytical level. 

It is worth noticing that we are introducing two new features in this fashion. In the original DMRG 
method of White the target state selected comprises just a few sites of the chain while now we are 
using the whole chain. On the contrary, in doing so we have to resort to variational or Fokker-Planck 
methods to handle the problem, while in the numerical DMRG the target state selected is exact for the 
particular size of the superblock chosen. 

4.1 Variational DMRG. 

In reference a new method based on the combination of perturbative and variational techniques 
was presented to study quantum lattice Hamiltonians. The general ideas of this method are illustrated 
in the example of the Ising model in an transverse field. The method relies on the choice of an 
exponential ansatz ?/>(A) exp[[/"(A)]?/>o> which is a sort of generalized lattice version of a Jastrow wave 
function. Perturbative and variational techniques are used to get successive approximations of the 
operator U(X). Perturbation theory is used to set up a variational method which in turn produces 
nonperturbative results. This method allows one to associate to the original quantum-mechanical 
problem a statistical-mechanical system defined in the same spatial dimension. In some instances these 
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statistical-mechanical systems turn out to be integrable, which allows one to obtain exact upper bounds 
to the energy. We shall briefly review this method hereby but for a detailed account we refer to 0. 

Let us suppose that we are given a hamiltonian of the form H(X) = Hq + A-Hi, where Hq has 
a nondegenerate ground state ipo and A is a coupling constant. In reference [J/J it was proposed to 
construct the ground state ip(A) of H(A) as 

oo 

VXA) = exp(^ X n U n ) *Po (4.83) 

n=l 

Solving perturbation theory in A to order v implies the knowledge of the collection of operators 
{C/ n } n= i i ... )l/ . Each operator U n consist in fact of a sum of "irreducible" operators Vj, 

U n = Y,PnjVi (4.84) 

Hence inserting ( 4783) into ( 4.83j ) and interchanging the order of the sums one arrives to a "dual" 
description of the ground state 

V>(A) = exp(]T aj(A)Vr) Vo (4.85) 
I 

where a/(A) = Y, n ^ n Pn,i- 

This expression suggests an alternative approximation to the ground state ip(A) which consists in 
choosing only a class of irreducible operators Vj whose weights oti are determined variationally. This 
was precisely the approach applied in Q to the Ising model in a transverse field. 

According to the perturbative-variational (PV) method @ we must first determine the form of 



the set operators {U n } by inserting the exponential ansatz (4.83) in the Schrodinger equation for 
H. Then these operators serve us to construct variational wave functions by inserting them back in 
the exponential ansatz. The perturbative equations that determine the U n for the lowest order in 
perturbation theory obey the equations, 

order A: ({H , U x ] + #i)Vo = £ (1) V>o (4.86) 

order A 2 : ([#„, U 2 ] + [H u f/i] + 1/2([[H , U x ], U x ]))^ = (4.87) 

where and E^> are the perturbative energies to first and second order. 
A solution of equations ( f4.86| ) and ( [4.87 ) in the case where Hq = — Y^j &j and H\ = — J2j a j a 
the following, 



^^E^K+i ( 4 - 88 ) 

5 

^ = ^{E a j a j+2 + E ^(^l+l + CT f-i)i (4-89) 

3 j 

In the following we shall make use of the simplest exponential ansatz based on the U\ operator only. 

Now let us start searching for the block density matrix using the variational method. Recall that 
the density matrix in Eq. ( f4.82| ) can be rewritten as a scalar product of two block states \iP)b an d 
\iP')b' given by their defining system-universe decompositon of the previous section: 

l^ = E E WW) ( 4 - 9 °) 

■ieB jeB c 
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W)B< = E E ^i'f\i')\f) ( 4 -91) 
i'eB' j/e B lc 

where i, i' are system- indexes, j, j' are universe-indexes, B, B' are system-blocks and B c , B' c are 
universe-blocks. The density matrix p is given by the scalar product of these two block states: 

P =B (4-92) 
for upon substitution of Eqs. flOOl) , ( ^9l|) we arrive at 



P = EE^IO#y#i'K')li / ) 

ij i'j' 

= EK)(%i' (4-93) 

with pjj/ as in Eq.( [4.82| ). 

Now let us consider the following variational state ansatz for the ground state |"0o) 1,2,3,..., N of the 
ITF model in a lattice with N sites: 

N 

|^o)l,2,3,...,iV = exp(^ a j a j+l) l°)l,2,3,...,iV (4.94) 
1 0=1 

where a is a variational parameter which is determined by minimizing the vacuum expectation value of 
the ITF Hamiltonian in this state, thereby making the parameter to become a function a = a(J/2T) of 
the coupling constant of the model; and 1 0) 1,2,3,..., = |0)i (g> . . . <g> \0)n with |0) being the ground state 
of the a x matrix. This exponential-variational ansatz constitutes part of a method called Perturbative- 
Variational method (PV) developed in for spin systems and in for fermionic systems. The PV 
method is essentially a cluster method which combines perturbative and variational techniques. Using 



Eq. ( 4.92 ) we construct the block density matrix out of this target state for a block containing 2 sites 
as: 

P PV = ^r-, \M*Wlv? (4-95) 
Z N (a) 

where Z^{a) is a normalization factor being the norm of the state which can be interpreted as the 
partition function of a certain associated statistical model [0]. As a matter of fact, in the large N limit, 
it turns out to be simply: 

Z N (a) = (ip \ipo) = (cosha)^, N > 1 (4.96) 

Inserting the ansatz (4.94) into fl4.95| ) we can express the density matrix in a more appealing form, 
namely, 

P PV = tAt-, i,2(0| exp ^(afa z 2 + a'\a'\) zf_ 2 {a, h, ti)\0) Vj2 , (4.97) 
Z N {a) 2 

where we have defined the following quantities, 



Z$l 2 (a,/ l ,/ l O : =3,..^(0|exp (a & 3 ^-l +1 +^+^)| 0)3i ^ jiV (4 . 98) 



h:= 2(^ + ^2) (4.99) 
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n 



(4.100) 



For the time being, it is convenient to shift N — 2 — > N in order to make expressions easier (at the end 

we shall come back to the correct value). 

It is possible to recast zff_ 2 (a, h, ti) into the following form, 



iV-l 



Z^>(a, h, ti) = 1 ,... )JV <0| exp(a ]T a z a z +1 + ha\ + tia z N )\0) h ... 



N 



N-l 



(4.101) 



{o-lv.Cjv} j = l 

Now we can recognize this equation as the partition function for the Ising model on an open chain of N 
sites subject to an external magnetic field applied only at the ends of the chain. This partition function 
can be worked out exactly using standard transfer matrix calculations yielding the result: 



,(o) 



1 



{[cosh(/i + h!) + cosh(/i — h')] (2 cosh a 



,N-1 



N 2^ 

+ [cosh(/i + ti) - cosh(/i - /i , )](2sinha) Af " 1 } 
In the N — > oo limit in which we are interested in, it further simplifies to: 

1 



7(0) 

J N^oo 



(cosh a) [cosh(/i + ti) + cosh(/i — h')] 



(4.102) 



(4.103) 



Inserting now Eq.(4.103) in Eq.( 4.97| ) we arrive at the following expression for the PV block density 
matrix: 



1 



o 



2(cosh a) 3 1 



<0|exp-(<7f<7 2 + aV2) 



x [cosh -{pi + a z 2 + a'\ + a%) + cosh -{a\ - a z 2 + a'\ - a' z 2 )]\0) V)2 ' 



(4.104) 



This is a nice result. Observe that the piece i, 2 (0| exp i^i^ exp^ 1<T |0) 1^,2' corresponds to a density 
matrix of a pure state p = \4>)y t 2' i,2(<t>\ ls the projection of the target state tpo onto the block (1,2). 
The extra terms in Eq. ( 4.104 ) are the novel features that the DMRG(PV) method brings about. 
To proceed further and give p PV a simple matricial form, it is convenient to change basis from eigenstates 
|0), |1) of a x to eigenstates |+), |— ) of a z . The notation is, 



In the new basis {|+) 



1+) 



1 



ID 



1 

V2 



-)} the components of p PV are: 

1 



exp — (crier? + 10" 2) 
a\a' 2 a 1 a 2 8(cosha) 3 P 2 v " " 



x [cosh ^((Ji + fj 2 + a\ + cr' 2 ) + cosh ^(01 - <7 2 + cr'i - cr' 2 )] 



(4.105) 
(4.106) 



(4.107) 
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Now p PV takes the following matricial form in the basis {++, 



-,+- 



" + }: 



/ e Q (cosha) 



„PV 



1 



4(cosh a) 



e a (cosh a)" 



cosh a 
cosh a 



cosh a 
V cosh a 

We can readly check that p PV is normalized: 



cosh a \ 
cosh a 

cosh a e _a (cosha) 2 e~ a 

cosh a e~ a e _Q (cosha) 2 J 



trp PV = 1 (4.109) 

Next step in the DMRG algorithm is to diagonalize the density matrix ( |4,108 ) and truncate to the 
largest ones; in this case the truncation is to two states to be denoted by \A) and \B) and will play the 
parallel role of the states \G) and \E) for the block Hamiltonian Hb of Section 3. 
We do not need to make a "blind" diagonalization of this 4x4 matrix for we may take advantage of 
what we have learnt in Section 3 about the eigenstates of the Hb in the ITF model. Then, the largest 
eigenvector say \A) will be in the even sector {|00}, |H}}, while the next to the largest one, say \B), 
will be in the odd sector {|01), |10)}. According to this analysis, we may write those states as, 



(4.108) 



\A) = x o|00) + a;ii|ll) 



(4.110) 



|B)=x 01 |01)+xio|10) (4.111) 
where xqq, x\\ , xqi , x%q are the components to be determined. Expressing the states 1 00) , |11) ... in the 



basis of {|+), |— )}, the diagonalization of the density matrix ( 4.108|) yields the following eigenvalues: 

1 



Wl 



4(cosh a) 3 
1 

4(cosh a) 3 



cosha(l + cosh a) + y cosh a(l + cosh a) 2 — sinh a] 



cosha(l + cosh^ a) — \J cosh^ a(l + cosh 2 a) 2 — sinh 4 a] 



1 



w 3 



W4 



4(cosh a) 3 
1 



e a sinh 2 a 



e a sinh 2 a 



4(cosh a) 3 

For a small it is easy to see that the eigenvalues are sorted according to 

U)Q > wz > w± > W4 



(4.112) 
(4.113) 
(4.114) 
(4.115) 

(4.116) 



For arbitrary values of a, which in turn amounts to arbitrary values of the coupling constant g = J/2T 
due to the variational equations to be given bellow, we have plotted these 4 eigenvalues in Fig. 6 
observing that there are not level crossings in the whole range of variation of a and that the sorting in 
Eq. fl4.116D holds all over, not just for small a. It is very important for our DMRG(PV) method to work 
properly that this property holds up, for when we truncate to the eigenstates of the largest eigenvalues, 
wo and ws, the physics will not change qualitatively when varying the coupling constant g. 
Moreover, it is also possible to show that the eigenvectors \A) (of w\) and \B) (of w^) are given by: 



|00) +a py |ll 



l + (a 



PV\2 



(4.117) 



17 



IB) = « (4.118) 
with the function a PV now the following function of the parameter a (after some tedious algebra), 



, M y' cosh a(l + cosh a) 2 — sinh a — 2 cosh a 
sinha(l + cosh 2 a) 



a = : : 2—s ( 4 ' 119 ) 



Notice that these states which now come from a DMRG(PV) analysis have the same form as the states 
\G), \E) of the block Hamiltonian Hg Eqs.( |3.37| )-( |3.38| ), the difference being in the dependence of the 
function a upon the coupling constant. This means that we have again the same structure as in the 
BRG analysis where the intertwiner operator Tq was fully determined, and consequently the whole RG 
procedure, by a single function a = a(g) of the coupling constant. To obtain a PV = a PV (g) we need 
the variational equation relating g with a. This can be found in |?j and is given by, 

<7=^ = tanha (4.120) 

Observe that when J/T <C 1 then a = J/2T and thus a PV ~ a/2 = J /AT. Now it is possible to 
eliminate the intermediate parameter a between Eqs. ( [4.119 ) and ( 4. 120| ) yielding the desired formula, 



'1- g 2 + g -r- (1- g 2 ) 

a PV ( 9 ) = 1 - 4 ' — (4.121) 

5(1 - \) 



Therefore we may appreciate that this function shares the same qualitatives properties as clbrg Q3.42 ) 



does, Eq.( p.56 ) , namely it goes to zero linearly when g — ► and it is bounded below 1, i.e., < a PV < 1 



for < g < 1. It addition, it has a singularity at g = \[2 which prevents the extension of this method 
beyond that singular point. The origin of this singularity is due to the variational nature of the method 
(see j7j). Nevertheless the critical region lies within the region of applicability of the PV method. 
Furthermore, it is also possible to define a variational DMRG method valid for the whole range of 
variation of the coupling constant. To do this, we simply recall that for small a the coupling constant 
depends linearly on the variational parameter, 

g ~ a, for a small (4.122) 
Thus, we may define another function say a (<?) by simply substituting a by g in a(a ) (pd9|) , 



/M , y'cosh g(l + cosh g) 2 — sinh g — 2coshg 
sinh g(l +cosh 2 g) 



a^ v =-S : g-r 1 ~ ( 4 - 123 ) 



In this fashion, a PV ' (g) shares the same properties with cibrGi without any singularity in the range of 
g. In fact, a PV (g) has a horizontal asimptota at 1 as qbrg does (see Fig 7). 



4.2 Fokker-Planck DMRG. 

This is another approximated version for preparing the target state to be projected onto the block- 
system in order to construct another analytical DMRG method based upon a blocking procedure. The 
details of how the Fokker-Planck (FP) method is applied to construct an approximate version of the real 
ground state of the ITF model are given in the original paper ||]. Briefly stated, what the FP method 
does is to start with an exponential ansatz as in the variational method Eq. ( |4.94j ), but instead the 
parameter a is fixed by demanding that a certain Fokker-Planck Hamiltonian H FP (to be determined 
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along the way) satisfies the eigenvalue equation to a certain order v in the perturbative expansion of 
the parameter a, that is, 



H FP (a)^ FP (a) = E FP (a)* FP (a) (4.124) 

such that the exact ITF Hamiltonian Hjtf an d its Fokker-Planck approximation version H FP ( a) 
differ in operators say Vi which involve interactions between lattice sites a v + 1 distant apart or larger. 
Squematically, 

H ITF - H FP {a) = J2 C/(a)V> (4.125) 

I>u 

Correspondingly, the FP-energy E FP (a), although incorporating non-perturbative effects, should agree 
with the exact ground state energy up to order v + 1 in a. In a sense, this gives the best "exact" 
aproximation to the Hamiltonian Hjtf to order v in perturbation theory |l. 

It is possible to show that these conditions fix the relationship between a and g. To lowest order this 
is given by [§, 

1 

c/ = -sinh(2a) (4.126) 

Again we see that for small coupling constant a and g are equal as in the variational method. 

In order to obtain the Fokker-Planck function, say a FP (g), we first notice that as in this FP method 
we start with the same exponential ansatz Q4.94j ) as in the variational method, the same function a(a) 
in Eq. ( [4,1 19j ) is valid here. The new feature is that we have to use the relation ( f4,126j ) now to express 
the function a FP (g) in terms of the coupling constant. After some tedious algebra we arrive at the 
following expression: 



i j> v ,. + 26 5 2 + 4^) + (8 + 6( 7 2 )v / TTV-2(l + yTTV) 

g{3 + a/1 + 4# 2 ) 

This function contains all the information which upon inserted in the intertwiner operator To gives rise 
to what we denote by a DMRG(FP) method. By looking at Fig. 7 we notice again that a FP (g) has the 
same qualitative properties as cibrg does. 

In this section we have introduced 3 functions a(g) namely, a PV (g), a PV (g) and a FP (g), related 
to different analytical realizations of the density matrix RG ideas. It is our purpose now to check the 
goodness of those methods by comparing their predictions for the critical exponents with the exact 
values already found by the standard BRG method in Sect. 3. 

In Fig. 7 we have plotted the 4 functions a{g) along with the universal function a(g c ) introduced in 
Eq. ( 3,58; ) whose cuts with the functions a(g) gives the predictions on the location of the critical point 



g c for each method. Recalling that the exact value is g^ xact = 1/2, we see that the closest value to this 
one is produced by the Fokker-Planck version of the DMRG method, even better than the standard 
BRG. Nevertheless, we may appreciate from Fig. 7 that the 4 methods lie rather close to one another 
within the critical region, the major differences being present off criticality when entering the strong 
coupling region. The particular values of g c are gathered in Table 2. 

Another interesting function to be plotted is the beta function f3(g) obtained for each method 
according to the analysis of Sect. 3, Eq. (3.61). We show these results in Fig. 8 where we observe that 



the 3 new beta functions introduced in this section by means of variational and Fokker-Planck DMRG 
methods have the same qualitative behaviour as the standard BRG beta function of Sec. 3. We know 
that in particular this means that the unstable character of the fixed point g c is preserved by these new 
methods. Moreover, we notice again that in the critical region the differences are small, namely, the 
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cut with the g-axis and the slope at the corresponding g c . These two latter properties are related to 
critical exponents. 

As far as the critical exponents is concerned, we have collected in Table 2 the exponents computed 
previously for the standard BRG method: correlation length exponent v, dynamical exponent z, mag- 
netic exponent f3 and the gap exponent s. Notice first that the 3 new DMRG methods also satisfy the 
scaling relation, 



which we know does not hold for the exact solution of the ITF model. 
From Table 2 we arrive at the following conclusions, 

• The best correlation length exponent v is provided by the DMRG Perturbative- Variational 
method. 

• The best dynamical exponent z is provided by the DMRG Fokker-Planck method. 

• The best magnetic exponent f3 is provided by the DMRG Perturbative- Variational method. 

• The best gap exponent s is provided by the DMRG Fokker-Planck method. 

From these results we may draw the conclusion that the RG methods based on block density matrix, 
either variational or Fokker-Planck, provide an improvement respect to the standard BRG methods, 
though it is not a major improvement. One of the reasons why this improvement is not as good as 
the numerical results obtained by White [flj relies on the fact that we have just kept 2 states in our 
analytical DMRG method while in the numerical treatment quite a lot states are kept. 



5 Conclusions 

We have presented in this paper an analytic formulation of the recently proposed Density Matrix RG 
method ||] . This method was originally developed in a numerical fashion mainly because it relies on the 
Wilsonian procedure of enlarging lattice sizes in the real-space RG. As this Density Matrix RG method 
has become a powerful tool to compute static and dynamic properties of quantum lattice systems at 
zero temperature, we find interesting to devise an analytic formulation of this method to be tested 



against the old procedures based upon the standard Block RG method introduced by Drell et al. [1C] 
to study QCD. 

The new feature of the DMRG method is its ability to take into account the unavoidable interaction 
between the block selected for truncating states and the rest of the lattice. This feature is the more 
relevant the more strongly correlated is the system under study, such as lattice QCD and strongly 
correlated electrons in High-T c materials. 

We have been able to devise an analytic formulation of the DMRG method by combining the idea of 
"interacting blocks" (that is, "non- isolated") with two other approaches. One approach is the old idea of 
blocking "a la Kadanoff" instead of the Wilsonian procedure. To this end we have presented an unified 
formulation of the old BRG method (sect. 3) in terms of what we call the intertwiner operator. This 
facilitates the task of bringing together the Kadanoff blocking with the DMRG. The second approach 
has been to devise an analytical method to produce "target states" which are the basic ingredients in 
the DMRG calculations. To this purpose we have used two recently proposed methods to deal with 
quantum lattice Hamiltonians: the Perturbative- Variational method @ and the Fokker-Planck method 
[||]. We have coined the names DMRG(PV) and DMRG(FP), respectively, for the density matrix RG 
methods coming out of these two approaches 

2 After we finished this work we have been aware of references where those authors also propose new extensions of 
the DMRG method. 
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In order to facilitate the task of testing these two new DMRG methods, we have chosen the simple 
ITF model where previous work with the BRG method is available. The results of computing several 
critical exponents are described in Sect. 4 where we have seen that either of the new DMRG methods 
perform better than the old BRG, although the improvement we get is not as good as White's numerical 
DMRG for we keep a much lower number of states during the truncation among other reasons. 

In addition, there is another line of work in order to improve our analytical DMRG methods. This 
is via the exact diagonalization of the block density matrix corresponding to 4 lattice sites (starting 
from the exact ground state for 4 lattice sites) and then to proceed in the usual fashion to construct 
the corresponding new a(g) function which would allow us to carry the blocking method. We leave this 
posibility open for future work. 

Another reason to seek analytical formulations of the DMRG method is the possibility of general- 
izing the current one-dimensional algorithms to higher dimensions. We consider our DMRG(PV) and 
DMRG(FP) methods as a first attempt at this goal. As a matter of fact, it is possible to use them to 
prepare 2-dimensional target states in the ITF model and to proceed with a DMRG analysis. Moreover, 
as far as fermion systems is concerned, it is also possible to apply the new density matrix methods using 
the perturbative- variational techniques described in fl~8|| . 

Finally, the renormalization group method is one of the basic concepts in several branches of Physics 
and we believe that we are currently facing a reconsideration of the old renormalization group ideas 
which will certainly have implications in areas such as Field Theory, Statistical Mechanics and Con- 
densed Matter. 
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Table captions 

Table 1 : Block Renormalization Group Method BRG1 (ne,A/). 

Table 2 : Critical exponents for the ITF model according to different RG methods and exact 
solution. 
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Figure captions 



Figure 1: Pictorical decomposition of the Hamiltonian H into single-site part Hg and two- nearest- 
neighbour-site part Hss- 

Figure 2 : Block decomposition of the open chain into blocks with ub = 3 sites. 

Figure 3 : Pictorical representation of the block Hamiltonian Hb and the interblock Hamiltonian 
H BB for the ITF model. 

Figure 4 : Pictorical representation of the truncation procedure for the block and interblock 
Hamiltonians in the ITF model. 

Figure 5 : Lattice decomposition into "system"- and "universe" -parts. 

Figure 6 : The 4 eigenvalues wq, w±, W3 and W4 of the block density matrix corresponding to the 
variational DMRG method. 

Figure 7 : Plot of the functions a(g) according to the methods: BRG (solid line), DMRG(PV) 
(grey line), DMRG(FP) (dashed line), DMRG(PV'). 

Figure 8 : Plot of the beta functions (3(g) according to the methods: BRG (solid line), DMRG(PV) 
(grey line), DMRG(FP) (dashed line), DMRG(PV'). 
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Table 1: Block Renormalization Group Method BRG1 (riB,k') 



Steps of the BRG1 Method 

1) Blocking Transformation: H = Hb + Hbb 

2) Diagonalization of Hb 

3) Truncation within each block: To 

4) Renormalization of Hb and Hbb 
H' s , = TqHbT 

H 's>s> = T^HbbTq 

5) Iteration: Repeat 1) -> 4) for H' = H' s , + H' s , 



Table 2: Critical exponents for the ITF model according to different RG methods and exact solution. 
Method g c v z (5 s 



BRG 0.3916 1.4820 0.5515 0.4086 0.8173 

DMRG(PV) 0.3790 1.4073 0.5353 0.3767 0.7534 

DMRG(FP) 0.4011 1.5177 0.5647 0.4285 0.8570 

DMRG(PV') 0.3874 1.4579 0.5459 0.3980 0.7958 

Exact Solution 0.5 1 1 0.125 1 
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